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ABSTRACT 

We present multiepoch spectra of 13 high-redshift Type la supernovae (SNe la) drawn from the 
literature, the ESSENCE and SNLS projects, and our own separate dedicated program on the ESO 
Very Large Telescope. We use the Supernova Identification (SNID) code of Blondin & Tonry to 
determine the spectral ages in the supernova rest frame. Comparison with the observed elapsed time 
yields an apparent aging rate consistent with the 1/(1 + z) factor (where z is the redshift) expected 
in a homogeneous, isotropic, expanding universe. These measurements thus confirm the expansion 
hypothesis, while unambiguously excluding models that predict no time dilation, such as Zwicky's 
"tired light" hypothesis. We also test for power-law dependencies of the aging rate on redshift. The 
best-fit exponent for these models is consistent with the expected 1/(1 -I- z) factor. 

Subject headings: cosmology: miscellaneous — supernovae: general 
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1. INTRODUCTION 

The redshift, z, is a, fundamental observational quan- 
tity in Friedman-Lemaitre-Robertson-Walker (FLRW) 
models of the universe. It relates the frequency of light 
emitted from a distant source to that detected by a local 
observer by a factor 1/(1 -I- z). One important conse- 
quence is that the observed rate of any time variation 
in the intensity of emi tted radiation w ill also be propor- 
tional to 1/(1 + z) (see lWeinberill972l and Appendix[J|- 

Due to their large luminosities (several billion times 
that of the sun) and variability on short timescales 
(^ 20 d from explosion to peak luminosity; iRiess et al.l 
119991 : IConlev et al.l[200i l. Type la supernovae (SNe la) 
are ideally suited to probe these time-dilation effects 
across a large fraction of the observable universe. The 
suggestion to use supernovae as "cosmic clocks" was pro- 
posed by Wilson more than six decades ago (Wilso^ 
|1939( ) and teste d on light c urves of low-redshift SNe la in 
the mid-1970s (|Rustlll974[ ). but only since the mid-1990s 
has this effect been unambiguou sly detected in the light 
curves of high-redshi ft objects (iLeibundgut et al.] 119961 : 
iGoldhaber etaII[200lD . 

These latter studies show that the light curves of dis- 
tant SNe la are consistent with those of nearby SNe la 
whose time axis is dilated by a factor (1 -I- z). However, 
there exists an intrinsic variation in the width of SN la 
light curves which is related to their peak luminosities 
(Phillips 1993) . such that more luminous SNe la have 
broader light curves (Fig. [T]). This "width- luminosity" 
relation is derived using low-redshift SNe la whe re the 
time-dilation effect 



if any, is negligible (f hiUips 
Hamuv et al.'1995 HRiess et al.lll995i|Phillips et all 



1993 



19991 : 



Jha et al. 2007). 



It is problematic to disentangle this intrinsic variation 
of light-curve width with luminosity and the effect of time 
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Fi g. 1. — B oloinctric light curves of 5 low-redshift SNe la taken 
from lStritzm ger (2005) (from top to bottom: SNe 1991T, 1999ee, 
1994D, 1992A, 1993H, and 1991bg). More luminous SNe la have 
broader light curves. SN 1991bg is an example of intrinsically 
subluminous SNe la (maximum Lbol < lO^Lsun), which are less 
likely to be found at high redshifts. 

dilation. To directly test the time-dilation hypothesis one 
needs to accurately know the distribution of light-curve 
widths at z ~ and its potential evolution with redshift, 
wh ether due to a selection effect (not taken into account 
by iGoldhaber et al.| |2001[ ) or an evolution of the mean 
properties of the SN la sainple wi th redshift — as possibly 
observed bv 'Howel l et al.l (|2007( ). Moreover, one needs 
to probe sufficiently high redshifts (z ^. 0-4, a s done by 
iLeibundgut et al.lll996HGoldhaber et al.ll2001h such that 
the observed widths of the SN la light curves are broader 
than the intrinsic width of any nearby counterpart. 

Furthermore, one might argue that at high redshift we 
are preferentially finding the brighter events (akin to a 
Malniquist bias). Such a selection effect would produce 
a spurious relation in which there would be broader light 
curves at higher redshifts, without any time dilation. 

The spectra of SNe la provide an alternative and a 
more reliable way to measure the apparent aging rate of 
distant objects. Indeed, the spectra of SNe la are remark- 
ably homogeneous at a given age, such that the age of a 
SN la can be determined from a single spectrum with an 
accuracy of 1-3 d ; — with no reference to its co rrespond- 
ing light curve iRJess et all Il997t iHowell et al. 2005; 
iBlondin fc Tonrvll2007f l. More importantly, the spectra 
of SNe la spanning a range of luminosities (and hence 
different in trinsic light-curve wid ths) evolve uniformly 
over time (jMatheson et al.l l2008f ) . The use of spectra 
thus avoids the degeneracy between intrinsic light-curve 
width and time-dilation effects. While there are some 
notable examples of deviatio ns from homoge neity in sev- 
eral SNe la (e .g., SN 200 0cx.lLi et al.l[2001al: SN 2002cx, 



Li et al.ll2003l: SN 2002ic. lHamuv et a l!''200g: SN 20G3fg, 
Howell et al.ll2006t SN 2006gz,|Hicto et al. 2007), these 



outliers arc readily identifiable spectroscopically through 
comparison with a large database of supe rnova spectra 
(see Section [2] and IBlondin fc Tonrvl [20071 ). 



As of today there are two published examples of ag- 
ing rate measurements usi ng spectra of a s ingle SN la 
(SN 1996bj a t z = 0.574 iRiess et al.l ll997l: SN 1997ex 
at z = 0.362. iFolev et al.ll2005f) . In both cases, the null 
hypothesis of no time dilation is excluded with high sig- 
nificance (> 95%). 

In this paper we present data on 13 high-redshift 
(0.28 < z < 0.62) SNe la for which we have multiepoch 
spectra . We use the Supernov a Identification (SNID) 
code of IBlondin fc Tonrvl ()2007[ ) to infer the age of each 
spectrum in the supernova rest frame, from which we de- 
termine the apparent aging rate of each SN la. These ag- 
ing rate measurements are then used to test the 1/(1-1-2:) 
time-dilation hypothesis expected in an expanding uni- 
verse. The data enable us for the first time to directly 
test the time-dilation hypothesis over a large redshift 
range. 

This paper is organized as follows. In Section [2] we 
explain how one can determine the age of a SN la based 
on a single spectrum and present the SNID algorithm 
used for this purpose. The aging rate measurements are 
presented in Section [3l and the time-dilation hypothesis 
(amongst others) is tested against the data in Section JH 
with the help of model selection statistics (information 
criteria). Conclusions follow in Section [5l 

2. DETERMINING THE AGE OF A SN la SPECTRUM 

2.1. SN la Spectral Evolution 

The spectra of SNe la consist of blended spectral lines, 
with a profile shape characteristic of stellar outflows. 
This line profile (also known as a "P Cygni" profile) con- 
sists of an emission component symmetric about the line 
center, and an absorption component that is blueshifted 
by the ^ 10,000 km s~^ expan sion velocity of the SN 
ejecta (jPinto fc Eastman! |2000( ). The expansion also 
causes a Doppler broadening of both components, such 
that a typical spectroscopic feature in SN la spectra has a 
width of ~ 100 A. As the ejecta expand, the photosphere 
recedes in the comoving frame of the supernova, such 
that the spectra probe deeper layers of the ejecta with 
time. Given the homologous nature of the expansion (ve- 
locity proportional to r adius), and the chemical stratifi- 
catio i i in the SN ejecta (iNomoto et"al]|1984l : iStehle et al.l 
l2005t iMazzali et al.l 120081 ). deeper lavers correspond to 
lower expansion velocities and an increased abundance 
of iron-peak elements. The impact on the spectra is 
twofold. First, the blueshift of SN la spectral lines de- 
creases with time (by as much as ^ 1000 km s~^ per day; 
iBenetti et alll2005l: IBlondin et al.ll2006[) . Second, due to 
the varying chemical composition at the photosphere, the 
relative shapes and strengths of spectral features change 
on a timescale of days. 

This complex spectral evolution is nonetheless pre- 
dictable to a large extent. At a given age, the spectra 
are remarkably homo geneous among d ifferent "normal" 
SNe la. According to lLi et al] (|2001b[ ). these constitute 
~ 65% of the local SN la sample, the rest consisting 
of intrinsically subluminous (~ 15%) or overluminous 
(~ 20%) events, whose spectra show deviations from 
those of normal SNe la. Subluminous SNe la are less 
likely to be found at high redshifts; in fact, no such object 
has been spectroscopically confirm ed in any high-redshift 
supernova search t o this day fe.g.. lMatheson et al.ll2005l : 
I Howell et al.l l2005l) . In what follows we only consider 
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Fig. 2. — Standard (light gray) and maxiniuni (dark gray) devi- 
ation from the mean spectrum (black line) for the 22 low-redshift 
SNe la for which we report an aging rate measurement (see Sec- 
tion |4jl, at four different ages — given in days from B-band max- 
imum light. A low-order curve has been divided out from each 
spectrum to reveal the relative shapes and strengths of the various 
spectroscopic features. 

normal and overluminous SNe la. In a separate paper 
we will show that none of the SNe la in the high-redshift 
sample presented here (see Section |4]) has a spectrum or 
light curve consistent with the subluminous variety of 
SNe la. 

The spectroscopic homogeneity of SNe la holds even 
when we consider both normal and overluminous objects 
in a representative sample of nearby events. In Fig. [5] we 
show the mean spectrum for the 22 low-redshift SNe la 
for which we report an aging rate measurement (see Sec- 
tion m, at four different ages. While there is an in- 
trinsic spectral variance among these different SNe la 
— some spectroscopic features correlate with luminosity 
(e.g., [Nugent et al. 1995; Mathcson et al. 2008), the av- 
erage deviations from the mean spectrum are small, and 
all spectra evolve in a similar manner over the course of 
several days, independent of light-curve width. 

Both the homogeneity and rapid evolution of SN la 
spectra enables an accurate determination of the age of 
a single spectrum. We explain how this is achieved in 
practice in the following section. 

2.2. The SNID Algorithm 

Given a large database of finely time-sampled SN la 
spectral templates, we can determine the age of a given 
input spectrum by finding the best-match template(s) 
in the database. T here are several stand ard tech- 
niques to do this (see iBlondin fc Tonrvl '2007 for a re- 
view) . In this paper, we use an implementatio n of the 
correlation techniques of iTonrv fc Davia (|1979f ): SNID 



(jBlondin fc Tonrvl I2007D . SNID automatically deter- 
mines the type, redshift, and age of a supernova spec- 
trum. We refer the reader to that paper for a more de- 
tailed discussion. 

The redshift of the input spectrum is a free parameter 
in SNID, although it can be fixed to a specific value. 
Comparison of the SNID redshifts with those determined 
from narrow emission and absorption lines in the host- 
galaxy spectra (t ypically accurate to < 100 km s^-'^; see 
iFalco et al.l 1X9991) yields a dispersion about the one-to- 
one corr espondence of only a^ y 0.005 out to a redshift 
z « 0.8 (|Blondin fc Tonrvll200l . 

Similarly, comparison of the SNID ages with those de- 
termined using the corresponding light curves yields a 
typical accuracy < 3 d, comp arable to other algorithms 
(Riess et ai1ll997HHowell et al. 2005 ). However, the age 
error is systematically overestimated. In this paper, we 
estimate the error as follows: each spectrum in the SNID 
database is trimmed to match the rest-frame wavelength 
range of the input spectrum, and is correlated with all 
other spectra in the database (except those correspond- 
ing to the same supernova). The age error is then given 
by the mean variance of all template spectra whose SNID 
age is within one day of the initial estimate. 

The success of SNID and similar algorithms lies pri- 
marily in the completeness of the spectral database. In 
Fig.[3]we show the age distribution of the SN la templates 
used in SNID for this paper (these do not include sub- 
luminous SNe la). This database comprises 959 spectra 
of 79 low-redshift (z < 0.05) SNe la with ages between 
— 15 and +50 d from maximum light. The spectra are 
taken from the literature, from public databases (such as 
SUSPECT23 or the CfA Supernova Archive^"), or from 
a set of unpublished spectra from the CfA Supernova 
Program. A full reference to all spectra in the SNID 
database is given by Blondir Tfc Tonrvl (120071 ). It is im- 
portant to note that each template spectrum is shifted 
to zero redshift and that each template age is corrected 
for the expected (1 + z) time-dilation factor. Because all 
the template SNe la are at low redshift {z < 0.05), this 
is a very small correction, and we will see in Section [H 
that this has no impact on the aging rate measurements. 
Thus, SNID determines the age a supernova would have 
at z = — that is, in the supernova rest frame. 

The number of SNe la shown in Fig. [3] is large enough 
that we can select a subsample (shown as an open his- 
togram) on which to conduct age determinations and ag- 
ing rate measurements on low-redshift SNe la. The size 
of this subsample is set by the requirement that remov- 
ing it from the SNID database would leave a sufficient 
number of templates i n a given age bin for a reliable 
age determination (see Blo ndin fc Tonrvl '2007) . It was 
also chosen to include a sufficient number of intrinsically 
overluminous SNe la: indeed, there are five such SNe la 
(SNe 1997br, 1998ab, 1999dq, 1999gp, and 2001eh) in 
this subsample, accounting for ~ 20% by number of ob- 
jects and spectra. For these specific tests, the templates 
corresponding to the 22 SNe la in this subsample are 
temporarily removed from the SNID database to avoid 
biasing the age determination. 

We deliberately restrict this subsample to ages between 

^•^ http://bruford.nhn.ou.edu/~suspect/indexl.html 

^■^ http://www.cfa.harvard.edu/supernova/SNarchive.html 
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Fig. 3. — Age distribution of the 79 SN la templates used in 
SNID (hatched histogram). There are a total of 959 spectra with 
ages between —15 and +50 d from S-band maximum light. The 
open histogram shows the subsample of 145 spectra from 22 SNe la 
(restricted to ages between —10 and 30 d from maximum) for which 
we report an aging rate measurement (see Section |4j . 

— 10 and +30 d from maximum light. Before —10 d, the 
number of spectral templates in the SNID database drops 
rapidly, and the age determination is inaccurate. Past 
+30 d, the spectra of SNe la do not evolve as rapidly as 
around maximum light, and the age determination is less 
precise (jBlondin fc Tonrvll2007t) . 

3. AGING RATE AT HIGH REDSHIFTS 
3.1. Spectroscopic Data 

Our aging rate measurements at high redshifts are 
based on a sample of 35 spectra of 13 SNe la in 
the redshift range 0.28 < z < 0.62. These include 
previously published dat a by the High-Z Supernova 
Search Team (SN 1996bi. iRiess et al.l ll997D. the S uper- 
nova Cosmol ogy Project (SN 199 7ex, ]Folev e t al."2005l: 
SN 2001go, iLidman et all [2001 . and the ESSENCE 
proje ct (SNe 2002iz, b027, and 2003js, iMatheson et all 
|2005| ). For SN 2001go we present our own reductions of 
the three epochs of spectroscopic data obtained from the 
ESO Science Archive Facil ity^^, as only th e first spec- 
trum was published by Lid man et al.l ()2005l ). The spec- 
tra of SN 04D2an (the highest-redshift SN la in this sam- 
ple) were obtained by members of the Supernova Legacy 
Survey (SNLS), and will be published as part of a larger 
sample of SNLS data by Stephane Basa and cowork- 
ers. SN 2006tk will be published alongside the complete 
ESSENCE supernova dataset in the near future. 

The other five SNe la (SNe 2006mk, 2006sc, 2007tg, 
2007tt, 2007un) are ESSENCE targets which were ob- 
served spectroscopically through two dedicated "Target- 
of-Opportunity" programs at the ESO Very Large Tele- 
scope^^. 

Details on the instrumental setup and data reduction 
can be found in the aforementioned references. The rest 
of the data will be presented more thoroughly in a sepa- 
rate publication. All of the data are shown in Fig. H) 

We use SNID to determine the redshift of each spec- 
trum. The redshift of a given supernova (^sn) is then 



SN 
(1) 


zqal 
(2) 


2^8 N 
(3) 


^SNID 
(4) 


1996bj 


0.574 


0.581 (0.005) 


0.580 (0.008) 
0.582 (0.008) 


1997cx 


0.361 


0.362 (0.002) 


0.362 (0.005) 
0.361 (0.004) 
0.362 (0.004) 


2001go 


0.552 


0.552 (0.005) 


0.552 (0.008) 
0.556 (0.008) 
0.550 (0.009) 


2002iz 


0.427 


0.425 (0.004) 


0.422 (0.006) 
0.428 (0.006) 


b027 




0.315 (0.003) 


0.315 (0.006) 
0.315 (0.004) 


2003JS 


0.363 


0.361 (0.003) 


0.359 (0.004) 
0.363 (0.006) 


04D2an 


0.621 


0.614 (0.006) 


0.608 (0.007) 
0.625 (0.011) 


2006mk 


0.475 


0.477 (0.003) 


0.479 (0.005) 
0.478 (0.007) 
0.474 (0.008) 
0.476 (0.006) 


2006SC 


0.357 


0.356 (0.004) 


0.355 (0.007) 
0.357 (0.007) 
0.356 (0.006) 


2006tJ< 




0.312 (0.003) 


0.312 (0.006) 
0.310 (0.003) 
0.315 (0.006) 


2007tg 




0.502 (0.004) 


0.503 (0.009) 
0.503 (0.008) 
0.502 (0.007) 


2007tt 


0.374 


0.376 (0.004) 


0.367 (0.008) 
0.379 (0.007) 
0.377 (0.005) 


2007un 


0.283 


0.285 (0.004) 


0.287 (0.006) 
0.285 (0.007) 
0.285 (0.005) 



^^ http://www.eso.org/sci/archive/ . 

26 Programs 078.D-0383 and 080.D-0477; PI: Jesper Sollerman. 



Column headings: (1) SN name. (2) Galaxy 
redshift (the typical error is < 0.001). (3) SN 
redshift, quoted as the error- weighted mean of 
the individual redshifts for each epoch. (4) SNID 
redshift for each epoch, in order of increasing age. 

computed as the error-weighted mean of the SNID red- 
shifts (zsnid) for each of its spectra (Table [T|). For ten 
SNe la, we also have a redshift determination from the 
host galaxy (zqal)- Comparison with the supernova red- 
shift shows an excellent agreement between the two mea- 
surements (better than 1%). For the three remaining 
SNe la (SNe b027, 2006tk, and 2007tg), only the SNID 
redshift is available, but the different redshift measure- 
ments for individual spectra all agree to within la, and 
we are confident about their accuracy. In what follows, 
we will use the galaxy redshift when available for the age 
and aging rate measurements. Given the excellent agree- 
ment between zqal and zqn, this choice has negligible 
impact on our results. 

3.2. Accuracy of Relative Age Determination 

An accurate determination of the rate of aging involves 
accurate knowledge of age differences. In what follows 
we test how well SNID determines differential ages using 
the subsample of low-redshift SNe la presented in Sec- 
tion 12.21 While the determination of absolute ages has 
no impact on the main result of this paper, we discuss 
their accuracy in Appendix |B] 

We determine the rest-frame age (igpcc) of each of the 
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Fig. 4. — Multiepoch spectra of the 13 high-redshift SNe la used in this study, binned to 10 A (gray). The vertical offset between each 
spectrum is for clarity only, and does not reflect differences in flux density (i<\; erg s~^cm~'^A~^) between them. In each plot, the age 
of the supernova increases downwards, and the observed time (in days) fr om the first sp e ctrum is indicated. Overplotted in black is a 
smoothed version using the inverse- variance- weighted Gaussian algorithm of IBlondin et all 120061) . 



145 spectra in the low-redshift subsample of 22 SNe la. 
We then compute the absolute age difference (Aigpec) 
between each unique pair of spectra corresponding to a 
given supernova. This amounts to 631 pairs. This age 
difference is then compared with the absolute observer- 
frame age difference (Atobs) for each spectrum pair. 
Since z « for this subsample, Aigpoc can be directly 
compared to Aiobs with no correction for time dila- 
tion. Given the restriction to ages between —10 and 
-1-30 d from maximum light in the low-redshift subsam- 
ple, Aispoc (and hence Atobs) is at most 40 d. 

The results are displayed in the upper panel of Fig. [51 
There is good agreement between Atgpoc and Aiobs , with 



a dispersion of only 2.0 d about the one-to-one corre- 
spondence. For Aiobs ^ 30 d, however, SNID systemati- 
cally underestimates the age difference by ~ 1.5 d. This 
is more apparent in the plot of residuals in the middle 
panel, it is mainly due to a systematic overestimate of 
rest- frame ages ispoc ^ — 7 d from maximum light, due 
to the lack of spectral templat es in the SNID databas e 
with similar ages (see Fig.l3landl Blondin fc Tonrvll2007f) . 
The lower panel of Fig. [5] shows the absolute fractional 
age difference vs. Aiobs- The quantity |Afobs/Atspcc| 
is a direct measure of the accuracy we can achieve for 
the aging rate determination. As expected, the frac- 
tional age difference decreases with increasing age differ- 
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Fig. 5. — Upper panel: Rest-frame age difference (Aispcc) vs. 
observer-frame age difference (Ato^s) for each spectrum pair for a 
given supernova. There are 631 such pairs, with a dispersion of 
2.0 d about the one-to-one correspondence. Middle panel: Residu- 
als in the upper panel vs. Ato^s- Lower panel: Ratio of Atspcc to 
Aiobsi again vs. Atobs- For '^iobs > 6 d, the fractional difference 
is less than 20%. 

ence. For Aiobs > 6 d, this difference drops below 20%. 
Tlie liigh-redshift data presented in tlie previous section 
span a sufficient range of observer-frame age difference 
that the aging rate determination is accurate. Note that 
the systematic underestimate of the age difference for 
A^obs ^ 30 d results in a negligible fractional difference. 

3.3. Aging Rate Determination 

The rest-frame age of each high-redshift SN fa spec- 
trum (ispec) is determined as outlined in Section [2.21 In 
each case, we fix the redshift to that determined in the 
previous section. The results are displayed in Table [21 
along with the corresponding observed date of each spec- 
trum (tobs)- However, since the aging rate determination 
depends on age differences (see previous Section) , we also 
report the observer-frame and rest-frame age from the 
first spectrum, respectively denoted Atobs and Atgpoc in 
Table HI 

We can then trivially compute the aging rate for each 
supernova. This is simply done through a least-squares 
fit of a line to Aigpec versus Atobs • The slope of the line is 
a measure of the aging rate, which should equal 1/(1 + z) 
in an expanding universe (see Appendix |X|. Were there 
no time dilation, the aging rate would equal one. The 
results are displayed in Fig. [6l 

We note that comparing the inverse of the slo pe in 
Fig. [S] (denoted "age factor" by iFolev et all I2005D and 
(1 + z) leads to asymmetric errors. The errors on the 
"age factor" [= (1 + z)] become highly non-Gaussian 
when the uncertainties of the individual age measure- 
m ents are large (> 1 d, as is the case in this paper, and 
in iFolev et al. 1 120051 for SN 1997ex), whereas the errors 



TABLE 2 

Observer-frame and rest-frame age 

differences 



SN 

(1) 



^obs 

(2) 



(3) 



(4) 



^f-spec 

(5) 



1996bj 



367.99 -2.2 
378.04 3.1 



1997ex 



815.08 
839.96 
846.03 



-1.6 

17.4 
21.2 



2001go 



2021.70 7.9 
2027.69 9.8 
2059.17 31.2 



2002iz 



2586.95 
2614.57 



-0.5 
17.6 



b027 



2589.95 
2616.57 



-3.5 

18.4 



2003JS 



2942.46 
2966.71 



-4.9 
12.5 



04D2an 



3026.20 
3032.20 



-2.5 
0.9 



2006mk 



4031.71 
4040.72 
4051.77 
4063.78 



-6.2 

-0.6 

7.3 

18.5 



2006sc 



4063.58 
4076.65 
4084.68 



0.9 

9.8 

13.4 



2006tk 



4089.57 
4100.57 
4103.59 



-8.8 
0.3 
2.9 



2007tg 



4381.75 
4391.65 
4405.61 



-6.1 
-0.5 
10.0 



2007tt 



4415.81 
4430.65 
4443.58 



-5.0 

6.1 

14.9 



2007un 



4441.61 
4451.60 
4460.58 



3.2 
11.1 

17.7 



0.00 
10.05 



0.0 
5.3 



0.00 0.0 
24.88 19.0 
30.95 22.8 



0.00 0.0 

5.99 1.9 

37.47 23.3 



0.00 
27.62 



0.0 

18.1 



0.00 
26.62 



0.0 
21.9 



0.00 
24.25 



0.0 

17.4 



0.00 
6.00 



0.0 
3.4 



0.00 

9.01 

20.06 

32.07 



0.0 

5.6 

13.5 

24.7 



0.00 
13.07 
21.10 



0.0 

8.9 

12.5 



0.00 
11.00 
14.02 



0.0 
9.1 

11.7 



0.00 

9.90 

23.86 



0.0 

5.6 

16.1 



0.00 
14.84 
27.77 



0.0 
11.1 
19.9 



0.00 

9.99 

18.97 



0.0 

7.9 

14.5 



Column headings: (1) SN name. (2) Julian date 
{JD) minus 2,450,000 at midpoint of observation. (3) 
SN rest-frame age in days from maximum light, de- 
rived from the cross-correlation with spectral tem- 
plates using SNID. (4) Observer-frame days from first 
spectrum. (5) Rest-frame days from first spectrum. 

on the aging rate [= 1/(1 + z)] are always Gaussian. 
This is illustrated in Fig. [7] using a Monte Carlo simu- 
lat ion of the age mea surements for SN 1997ex presented 
bv lFolev et al.l (|2005D . Using the same errors on the in- 
dividual age measurements, the distribution of the slope 
measurements is highly non-Gaussian in (1 -I- z) space, 
while it is normally distributed in 1/(1 + z) space. 

The individual aging rate measurements presented here 
alone reject models that predict no time dilation at a high 
significance (up to ~ 6cr) , and all (except for SN 2006mk) 
arc within Ict of the expected 1/(1-1-2;) factor. In the next 
section we combine all aging rate measurements (includ- 
ing those for the low-redshift sample) to test each hy- 
pothesis more thoroughly. 

4. TESTING THE 1/(1 + z) TIME-DILATION HYPOTHESIS 

We have determined the aging rate for the subsample 
of 22 low-redshift SNe la. We combine these aging rates 
with those determined for the 13 high-redshift SNe la of 
our sample (see Table[3]) to test the l/(l-|-z) time-dilation 
hypothesis. As noted in Section 13 21 and Appendix JBl 
these measurements rely on a database of SN la spectra 
whose ages have already been corrected for the expected 
1/(1 + z) time-dilation factor. However, since all SNe la 
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Fig. 6. — Comparison of rest-frame (Atspoc) and observer-frame (Aiobs) time from the first spectrum, for each of the 13 high-rcdshift 
SNe la in our sample. The abscissa and ordinate ranges are both set to [—3, +40] d in all cases. The slope of the best-fit line (solid line) 
gives a measurement of the apparent aging rate of the supernova, which is compared to the expected 1/(1 + z) value. The dotted line in 
each plot corresponds to Atspcc = Atobg. 

the aging rate measurements. 

All aging rate measurements are shown in Fig. [51 The 
solid line shows the expected 1/(1 -I- z) time-dilation fac- 
tor, while th e dashed line represents the "tired light" 
hypothesis of lZwickvl ()1929f ). According to this hypothe- 
sis, photons lose energy as they interact with matter and 
other photons in a static universe. The energy loss is pro- 
portional to the distance from the source, and causes a 
redshift in spectra as in an expanding universe. However, 
this hypothesis does not predict a time-dilation effect, 
and so the aging rate should equal one for all redshifts. 

As expected, the measurement of a time-dilation effect 
is more obvious at larger redshift, and the precision im- 
proves as the number and time span of spectra for each 
supernova increases. This latter effect explains why the 
aging rate measurements for SN 1996bj (z — 0.574) and 
SN 04D2an [z — 0.621) have a large associated error 
— despite being the two highest-redshift SNe la in our 
sample, since only two spectra separated by ~ 10 (for 








(X - n-^la-^ ; X= slope 



Fig. 7. — Monte Carlo results illustrating the advantage of work- 
ing in 1/(1 + z) space for time-dilation measurements. Solid lines: 
recovered slope (in standard deviations from the mean, n) in (l-l-z) 
(left) and 1/(1 - I- z) sp ace (right), using the SN age errors reported 
by Fo lev et al.l I I2005I ). The distribution is highly non-Gaussian in 
the former case. 

in the SNID database are at redshifts z < 0.05, the cor- 
rection is small (^ 1 d) and has a negligible impact on 
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Fig. 8. — Apparent aging rate versus 1/(1 + z) for the 13 high-redshift (0.28 < z < 0.62) and 22 low-redshift {z < 0.04) SNe la in our 
sample. Overplotted are the expected 1/(1 + z) time dilation (solid line) and the best-fit 1/(1 + z)' model (with b = 0.97 ± 0.10; dotted 
line and gray area). The dashed line corresponds to no time dilation, as expected in the tired-light model — clearly inconsistent with the 
data. The inset shows a close-up view of the low-redshift sample. These data are summarized in Table [3] 



SN 1996bj) and - 6 (for SN 04D2an) observer-frame 
days are available. 

A simple x^ analysis is sufficient to confirm what the 
eye sees: the hypothesis of no time dilation is not a good 
fit to the data (x^ = 150.3 for 35 degrees of freedom; 
see Table m, with a goodness-of-fit of ~ 0% (defined as 
GoF= r(j//2,xV2)/r(i//2), where T{iy/2,x^/2) is the 
incomplete gamma function and v is the number of de- 
grees of freedom) — namely, a null probability of ob- 
taining data that are a worse fit to the model, assuming 
that the model is indeed correct. The expected 1/(1 + z) 
time-dilation factor, on the other hand, yields a good 
fit to the data (x^ = 27.0 for 35 degrees of freedom), 
with GoF= 83.2%, and is largely favored over the null 
hypothesis of no time dilation (Ax^ ~ 123). 

This result holds (and in fact improves) when we con- 
sider only the high-redshift sample (see Table S]) . This 
works because the z = end of the aging rate versus red- 
shift relation (Fig. [8]) is fixed to unity by theory, so the 
low-redshift sample is not needed to anchor the theoret- 
ical curve at z « (although it is still used to calibrate 
the ispcc measurement). The low-redshift data alone do 



not enable us to distinguish between the two hypothe- 
ses, since the impact of time dilation is small at such low 
redshifts. 

In Fig. [5] a different view of Fig. [5] shows the distribu- 
tions of the ratio between the aging rate and 1/(1 -I- z) 
for both the low-redshift (open histogram) and high- 
redshift (hatched histogram) samples. Both distributions 
are within ~ 20% of a unit ratio, again validating the 
hypothesis of time dilation over a large redshift range. 
The apparent bias to lower values of the ratio for the 
low-redshift sample is not statistically significant, as the 
mean error on the aging rate is of order one bin size 
(< 0.1; seeTableEl). 

In what follows we test whether the data favor a non- 
linear dependence of the aging rate on redshift, namely 



aging rate 



1 



(1 + z)' 



(1) 



where fe is a free parameter. While Eq. [T] satisfies the 
same zero point as the two previous hypotheses (aging 
rate equal to 1 at z = 0), no model actually predicts such 
a dependence of the aging rate on redshift. Nonetheless, 
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TABLE 3 

Aging rate measurements 

z 1/(1 + z) Aging rate 





Low rodshift {z < 0.04) 




1981B 


0.006 


0.994 


1.099 


0.071) 


1989B 


0.002 


0.998 


1.036 


0.079) 


1992A 


0.006 


0.994 


1.040 


0.077) 


1994D 


0.002 


0.998 


1.011 


0.152) 


1996X 


0.007 


0.993 


0.886 


0.116) 


1997br 


0.007 


0.993 


1.120 


0.149) 


1998V 


0.018 


0.983 


1.047 


0.166) 


1998ab 


0.027 


0.974 


0.987 


0.151) 


1998dm 


0.007 


0.993 


0.778 


0.119) 


1998eg 


0.025 


0.976 


0.883 


0.085) 


1999cl 


0.008 


0.992 


0.870 


0.122) 


1999dq 


0.014 


0.986 


0.928 


0.055) 


1999ej 


0.014 


0.986 


1.147 


0.159) 


1999gp 


0.027 


0.974 


0.839 


0.060) 


2000fa 


0.021 


0.979 


0.884 


0.076) 


2001V 


0.015 


0.985 


0.973 


0.072) 


2001eh 


0.037 


0.964 


1.148 


0.107) 


2001ep 


0.013 


0.987 


0.925 


0.086) 


2002ha 


0.014 


0.986 


0.973 


0.134) 


2003cg 


0.004 


0.996 


1.006 


0.131) 


2003du 


0.006 


0.994 


0.909 


0.091) 


20061f 


0.013 


0.987 


0.941 


0.110) 




High redshift (z 


>0.2) 




1996bj 


0.574 


0.635 


0.527 


0.369) 


1997ex 


0.361 


0.735 


0.745 


0.076) 


2001go 


0.552 


0.644 


0.652 


0.062) 


2002iz 


0.427 


0.701 


0.655 


0.089) 


b027 


0.315 


0.760 


0.823 


0.092) 


2003JS 


0.363 


0.734 


0.718 


0.082) 


04D2an 


0.621 


0.617 


0.567 


0.341) 


2006mk 


0.475 


0.678 


0.753 


0.060) 


2006SC 


0.357 


0.737 


0.619 


0.121) 


2006tk 


0.312 


0.762 


0.835 


0.181) 


2007tg 


0.502 


0.666 


0.687 


0.102) 


2007tt 


0.374 


0.728 


0.718 


0.108) 


2007un 


0.283 


0.779 


0.759 


0.135) 
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Fig. 9. — Ratio of the aging rate to 1/(1 + ^) for all SNe la in 
Fig. \E\ Both the low-redshift (open histogram) and high-redshift 
(hatched histogram) samples are shown. 

small deviations from the expected l/(l + z) factor would 
have profound implications for our assumption of FLRW 
cosmology. 

Again, we performed a least-squares fit to the entire 
sample, and also to the individual high- and low-redshift 
samples (see Table H]) . The data constrain the b expo- 
nent to 10% (Icr), and yield b = 0.97 ± 0.10 for the en- 



tire sample (Fig. [Sj dotted line and gray region) and 
b = 0.95 ± 0.10 for the high-redshift sample. As ex- 
pected, the low-redshift sample alone is insufficient to 
constrain the free parameter {b = 3.18 ± 1.28). Nonethe- 
less, the samples that include the high-redshift objects 
have a best-fit value for b that is consistent with 6=1, 
and thus with the expected l/(l-|-z) time-dilation factor. 
Since this model has an additional free parameter, it 
is instructive to use information criteria to compare it to 
the simple 1/(1 -|-z) prediction. These model comparison 
statistics favor models that yield a good f it to the data 
with fewer parameters. As in lDavis et al.l (l2007f). w e use 
the Akaike information criterion (AlC: [Akaikelll97^ . For 
Gaussian errors (which is the case here, see Section [ 
this criterion can be expressed as 



AlC = x' 



2fc, 



(2) 



1.3 tei 



where k is the number of free parameters (jDavis et al.l 
|2007| ). Comparison of models simply involves comput- 
ing the difference in AlC (AAIC) with respect to the 
model with the lowest value for this criterion. A differ- 
ence in AlC of 2 is considered positive evidence against 
the model with the higher AlC , whereas a difference of 
6 is c onsidered strong evidence ()Liddlell2004l : iDavis etldl 
|2007| ). In the models considered here (see Table S]), the 
expected 1/(1 + z) time-dilation model has the lowest 
AIC (although this is not true for the low-redshift sam- 
ple), and we compute AIC differences with respect to 
that model. 

With AAIC = 1, we conclude that the information 
criteria do not provide positive evidence against a 1/(1-1- 
zf' dependence of the aging rate. The x^ per degree of 
freedom is also satisfactory for the samples that include 
the high-redshift objects. 

The two other models considered previously have no 
free parameters (A: = in Eq. [2]), hence AAIC = A^^, 
and the information criterion is reduced to a simple x^ 
test. 

5. CONCLUSION 

We have presented 35 spectra of 13 high-redshift 
(0.28 < z < 0.62) SNe la, which include previously un- 
pubhshed data from the ESSENCE and SNLS projects 
and from our own dedicated program at the ESO Very 
Large Telescope. Given the rapid and predictable evolu- 
tion of SN la spectral features with age, as well as the 
relative homogeneity of SN la spectra at a given age, one 
is able to determine the (rest-frame) ag e of a single spec- 
trum with a typical accuracy of 1-3 d (Riess et al.' 1997; 
^ley et al. 2005; Hook et all 120051 : [HowcU et al. 2005; 
Blondin & Tonrvl l2007[) . 

Using the Sup e rnova Identification (SNID) code of 
iBlondin &: Tonrvl (|2007l ). we determine the ages of each 
spectrum in the supernova rest frame. Comparison with 
the observed time difference between the spectra yields 
an apparent aging rate consistent with 1/(1 -I- z), as ex- 
pected in a homogeneous and isotropic expanding uni- 
verse. Moreover, the d ata unambigu ously rule out the 
"tired light" hypothesis (|Zwickvlll929[ ) in which photons 
lose energy as they interact with matter and other pho- 
tons in a static universe. 

The fact that the age determination is so accurate over 
a large redshift range also shows that the deviations be- 
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TABLE 4 
Time-dilation model comparison 



All SNo 



Model'' 



XVdof 



GoF 

(%) 



AAIC 



High-redshift SNe only 
X"/dof 



GoF 

(%) 



AAIC 



Low-redshift SNe only 
xVdof GoF AAIC 

(%) 



1/(1 + z) 

1/(1 + ^)" 
tired light 



27.0/35 

26.9/34 

150.3/35 



83.2 

80.2 

0.0 





1 

123 



3.6/13 

3.4/12 

123.4/13 



99.5 

99.2 

0.0 





1 

119 



23.4/22 
20.3/21 
26.9/22 



38.2 
50.0 
21.4 



The best-fit values for the b exponent in the second model are as follows. All SNe: 
0.97 ± 0.10; high-redshift SNe only: b = 0.95 ± 0.10; low-redshift SNe only: b = 3.18 ± 1.28 



tween spectra of low- and high-redshift SNe la in our 
sample are small. 

We also test for alternate dependencies of the aging 
rate on redshift, namely 1/(1 -I- z^, although these are 
not predicted by any model. Whether we consider the 
entire sample or only the high-redshift sample, the best- 
fit value for the b exponent is consistent with 5=1, and 
thus with the expected 1/(1 -t- z) factor. 

That these data provide a confirmation of the time- 
dilation factor expected in an expanding universe should 
be of no surprise. Nonetheless, previous u se of SN la light 
curves to test this hy pothesis (Lcibund gut et all 119961 : 
iGoldhaber et al.ll200lD are prone to the spread in intrin- 
sic light-curve widths and its possible variation with red- 
shift (which includes selection effects; see Section [1}. 

The data presented here are unique in that they en- 
able the most direct test of the 1/(1 -I- z) time-dilation 
hypothesis over a larger redshift range than has yet been 
performed. This hypothesis is favored beyond doubt over 
models that predict no time dilation. With more data, 
the focus will shift to testing more thoroughly the alter- 
native l/(l-|-z)'' dependence of the aging rate on redshift. 
Any significant deviation from 6 = 1 would have a pro- 
found impact on our assumption of a FLRW cosmology 



to describe the universal expansion. 
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APPENDIX 
A. TIME DILATION IN AN EXPANDING UNIVERSE 



In a homogeneous, isotropic expanding universe, the i nterval dr be tween two space-time events is given by the 
Robertson- Walker (RW) metric (lRobertsonlll935L Il936allbl : I Walked fl936l) . 



dr^ =c^dt^ -a^{t) 



dr'^ 



1 — kr"^ 



r^idO^ 



(Al) 



where c is the speed of light in vacuum, t is the cosmic time, {r,9,(j)) are the comoving spatial coordinates, k is the 
curvature parameter, and a{t) is the dimensionless scale factor. In what follows we assume the present-day value of 
the dimensionless scale factor uq — 1. 

Photons travel along null geodesies (dr^ =0). In what follows we consider radial null rays o nly ( d9 — d(f) — 0). For 
a photon emitted at time ti from an object located at (ri, 9i,(j)i) and observed at time to, Eq. |Al| implies 



*° cdt 



dr 



VI - fcr2 



fin)- 



(A2) 



Here we assume that the object from which the photon was emitted has constant coordinates (ri,0i,0i) such that 
/(ri), also known as the comoving distance, is time independent. Thus, for a photon emitted at time ti -|- Sti and 
observed at time to + Sto, Eq. |Alj also implies 



fin). 



(A3) 



'■*°+''*° cdt 

iti+sti ait) 

For small 6ti (and hence small 6to), the rate of change of the scale factor remains roughly constant and Eqs. |A2j and 
^3j imply 



6ta Sti 
ao a(^i) 



(A4) 
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Hence, a light signal emitted with frequency vi will reach us with frequency vq such that 

I/O Sti a{ti) 



vi Sto ao 



(A5) 



Using the standard definition of redshift, z ^ {Iq — h)/h = vijvii — 1, wc obtain a relationship between observed 
and rest-frame time intervals in a RW metric as a function of redshift z: 

A supernova at redshift z will thus appear to age (1 + z) times more slowly with respect to a local event at z « 0. 

The prediction of time dilation proportional to (l+z) is generic to expanding universe models, whether the underlying 
theory be general relativity (e.g., the Friedmann-Lemaitre- Robertson- Walker universe), special relativity (e.g., the 
Milne Universe), or Newtonian expansion. A point of confusion can occur in the special relativistic case for which the 
well-known time-dilation factor is given by 



7SR = 



-1/2 

(A7) 



which evidently differs from (1 + z). Thus it might be assumed that a special relativistic expansion can be distinguished 
from the FLRW universe using a time-dilation test^^. 

This is not the case. Special relativistic expansion of the universe assumes there is an inertial frame that extends 
to infinity (impossible in the non-empty general relativistic picture) and that the expansion involves objects moving 
through this inertial frame. The time-dilation factor from Eq. |A8| relates the proper time in the moving emitter's 
inertial frame ((5ti) to the proper time in the observer's inertial frame ((5io). To measure this time dilation the observer 
has to set up a set of synchronized clocks (each at rest in the observer's inertial frame) and take readings of the 
emitter's proper time as the emitter moves past each synchronized clock. The readings show that the emitter's clock 
is time dilated such that (5to = TSR^ii- 

We do not have this set of synchronized clocks at our disposal when we measure time dilation of supernovae in an 
expanding universe and therefore Eq. |A8| is not the time dilation we observe. We must also take into account an 
extra time-dilation factor that occurs because the distance to the emitter (and thus the distance light has to propagate 
to reach us) is increasing. In the time St^ the emitter moves a distance vit^ away from us. The total proper time we 
observe, (5to,tot, is (5io plus an extra factor describing how long light takes to traverse this extra distance (w(5io/c), 

<5io,tot = ^io(l + ^'/c). (A9) 

The relationship between proper time at the emitter and proper time at the observer is thus 

(AlO) 

(All) 

= (5ti(l-fz), (A12) 

which is identical to the GR time-dilation equation. 

Non-cosmological redshifts (i.e., not due to universal expansion) also cause a time-dilation effect described by 
Eq. |A6| . However, these additional effects from peculiar velocities and gravitational redshifts contribute random 
error only, and do not bias the measurements presented here. 

B. COMPARISON OF SPECTRAL AND LIGHT-CURVE AGES 

To test the accuracy of the age determination using SNID, we select the SNe la for which a well-sampled light curve 
is available around maximum light. Only the ESSENCE and SNLS SNe la in our sample have associated light curves 
for which we could dete r mine the date of maximum brightne s s (tma x). To do so we used the MLCS2k2 hght-curve 
fitting code of iJha et al.l ()2007t ). as done bv I Wood- Vasev et al.l (|2007D . This way we can determine the time difi'erence 
(in the observer frame) between maximum light (imax) and the time the spectrum was obtained (iobs)- We compare 
this time interval with the rest-frame age determined through cross-correlation with local SN la spectral templates 
using SNID (tspcc)- We expect a one-to-one correspondence between 

, ^obs ^inax /j-j -. n 

*LC = T-. (£>1) 

1 + z 
^"^ In fact, such an erroneous assumption was made by one of the current authors of lDavis &: Lineweaveri II2004I ) . 
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Fig. 10. — Upper panel: comparison of supernova rest-frame ages (in days from maximum light) obtained from cross-correlation with 
spectral templates (tspcc) and from fits to the light curve (4lc)- 145 age measurements for the subsample of 22 low-redshift SNe la are shown 
in gray. The dashed line represents the one-to-one correspondence between Jlc ^-iid ispcc- Middle panel: Age residuals, At = tspcc — tLC- 
We also indicate the standard deviation (cr) and mean residual (fi). Lower panel: Same as above, where each point has been corrected for 
the mean offset between tspoc and tLC for a given supernova. 

and tspec- 

The result is shown as black points in Fig.[Tni While the agreement is good, there is a mean systematic offset of —1.6 
d between igpcc and ^lCj as shown in the middle panel. If this offset were to affect only a subset of age measurements 
for a given supernova, the impact on the aging rate determination would be severe. To check this, we correct the 
spectral ages of a given supernova for the mean difference between ^lc and ispec- This "corrected" age residual, Atcorr 
is plotted in the lower panel of Fig.[TOl The mean residual drops to —0.1 d and the scatter decreases slightly. 

Since there are 2 to 4 ispec measurements for a given supernova, and only one measurement of imax, the source of the 
discrepancy between the spectral and light-curve ages is most likely due to the determ ination of the date of maximum 
using the light-curve fitter. Indeed, using a different light-curve fitter (SALT2; iGuv e t al. 2007) yields values for i^ax 
that differ from the MLCS2k2 measurements by more than one day in 9 out of 10 cases, and by more than two days 
for three objects (SNe 2003js, 2007tg, and 2007un). These discrepancies are due to a combination of differences in 
light-curve fitter algorithms and d ata quality (light-curv e sampling around maximum light and signal-to-noise ratio of 
each light-curve measurement; see iMiknaitis et aLll2007f) . 

Therefore, while there is a systematic offset between part of these different age determinations, this offset affects all 
measurements in a similar fashion and has no impact on the determination of the rate of aging. In fact, the main result 
of this paper (see Sectional) is completely independent of iLC; and hence of imax- Nonetheless, the comparison betwee n 
spectral and light-curve ages confirms the accuracy of age determination using spectra alone (jBlondin fc Tonrvll2007D . 
The age measurements for all the high-redshift SNe la in our sample are reported in Tabled 

In making the comparison we have implicitly assumed what we are trying to show, namely a time-dilation factor of 
(l + z). Accordingly, we also make the same comparison for our subsample of 22 low-redshift SNe la (0.002 < z < 0.04). 
At such low redshifts, the (1 -|- z) correction present in ^lc is negligible (the mean correction is ^ 0.06 d). The result is 
shown as gray points in Fig. [TOl The mean residual between ^lc and igpcc for this low-redshift sample is close to zero 
with a small scatter (a ~ 1.5 d), and unlike the high-redshift sample there is no significant systematic offset between 
the two age measurements. 

The age measurements presented in Table [5] also enable us to infer the date of maximum light for each supernova 
using spectra alone (corresponding to ispcc — 0). This way we are able to determine the time of maximum for the SNe la 
in our sample for which a well-sampled light curve was unavailable (SNe 1996bj, 1997ex, and 2001go; see Table [6]). We 
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TABLE 5 

Comparison of rest-frame light-curve and spectral 

AGES 



SN 
(1) 



tLC 
(2) 



''Spec 

(3) 



At 
(4) 



Ate 



(5) 



2002iz 


0.1 


1.1 


-0.5 


2.2 


-0.6 


2.4 


0.8 (1.5) 




19.4 


1.1 


17.6 


1.2 


-1.8 


1.6 


-0.4 (2.3) 


b027 


-2.4 


0.5 


-3.5 


1.8 


-1.1 


1.9 


-0.9 (1.7) 




17.9 


0.5 


18.4 


1.6 


0.5 


1.7 


0.7 (1.4) 


2003JS 


-3.2 


0.3 


-4.9 


1.6 


-1.7 


1.6 


0.3 (2.0) 




14.6 


0.3 


12.5 


1.2 


-2.1 


1.2 


-0.1 (2.3) 


04D2an 


-3.2 


0.9 


-2.5 


1.6 


0.7 


1.8 


0.2 (1.4) 




0.5 


0.9 


0.9 


1.3 


0.4 


1.6 


-0.1 (1.3) 


2006mk 


-3.6 


0.7 


-6.2 


1.0 


-2.6 


1.3 


-0.5 (2.7) 




2.6 


0.7 


-0.6 


2.2 


-3.2 


2.3 


-1.1 (3.3) 




10.0 


0.7 


7.3 


1.9 


-2.7 


2.0 


-0.6 (2.8) 




18.2 


0.7 


18.5 


1.8 


0.3 


1.9 


2.4 (0.9) 


2006SC 


0.9 


0.5 


0.9 


1.6 


0.0 


1.7 


0.9 (1.0) 




10.5 


0.5 


9.8 


1.4 


-0.7 


1.5 


0.2 (1.2) 




16.4 


0.5 


13.4 


2.2 


-3.0 


2.3 


-2.1 (3.2) 


2006tk 


-6.1 


0.5 


-8.8 


2.4 


-2.7 


2.5 


-0.8 (2.8) 




2.3 


0.5 


0.3 


2.0 


-2.0 


2.1 


-0.1 (2.2) 




4.6 


0.5 


2.9 


0.9 


-1.7 


1.0 


0.2 (1.9) 


2007tg 


-6.6 


0.9 


-6.1 


2.0 


0.5 


2.2 


0.2 (1.2) 




0.0 


0.9 


-0.5 


1.8 


-0.5 


2.0 


-0.8 (1.2) 




9.3 


0.9 


10.0 


1.5 


0.7 


1.8 


0.4 (1.3) 


2007tt 


-2.5 


0.6 


-5.0 


2.1 


-2.5 


2.2 


-0.1 (2.7) 




8.3 


0.6 


6.1 


1.4 


-2.2 


1.6 


0.2 (2.5) 




17.7 


0.6 


14.9 


2.2 


-2.8 


2.2 


-0.4 (3.0) 


2007un 


4.3 


0.4 


3.2 


2.2 


-1.1 


2.3 


0.1 (1.4) 




12.1 


0.4 


11.1 


1.3 


-1.0 


1.3 


0.2 (1.4) 




19.1 


0.4 


17.7 


1.4 


-1.4 


1.5 


-0.2 (1.7) 


Column 


headini, 


s: 


1) SN ne 


Lme 


(2) SN rest-frame age in 


days from 


maxim 


um 


light, der 


ivec 


from th 


c light curve. (3) 


SN rest-frame age 


in 


days from 


maximum light, 


derived from 


the cross-correlati 


on 1 


with spec 


tral 


template 


3 using SNID. (4) 



At: 



tsspcc and tLC- 



"^LC- (5) ^* corrected for the mean offset between 



TABLE 6 

Comparison of dates of maximum light 



SN 
(1) 



''max 

(2) 



Atn 



(3) 



(4) 



1996bj 




372.16 


(3.73) 




1997ex 




817.16 


(1.98) 




2001go 




2011.47 


(3.12) 




2002iz 


2586.83 (1.51) 


2587.71 


(3.21) 


0.88 (3.54) 


b027 


2593.09 (0.65) 


2594.20 


(1.90) 


1.11 (2.01) 


2003JS 


2946.80 (0.47) 


2949.29 


(1.67) 


2.49 (1.73) 


04D2an 


3031.36 (1.50) 


3030.61 


(1.82) 


-0.75 (2.36) 


2006mk 


4036.95 (0.96) 


4040.32 


(1.02) 


3.37 (1.40) 


2006SC 


4062.39 (0.64) 


4061.71 


(2.77) 


-0.68 (2.84) 


2006tk 


4097.56 (0.66) 


4100.13 


(0.98) 


2.57 (1.18) 


2007tg 


4391.62 (1.34) 


4391.34 


(1.58) 


-0.28 (2.08) 


2007tt 


4419.29 (0.83) 


4422.47 


(1.83) 


3.18 (2.01) 


2007un 


4436.12 (0.47) 


4437.15 


(3.12) 


1.03 (3.16) 



Column headings: (1) SN name. (2) JD - 2,450,000 of 
maximum light, derived from the light curve. (3) JD — 
2,450,000 of maximum light, derived from the spectra. (4) 

Aj. _ iSpcc ,LC 

i-it-max — ''max '•max. 

can also compare the dates of maximum as inferred from a fit to the light curve (inj^x) with those determined from 
the spectra alone (t^fx)- The results are also shown in Table [51 For four objects (SNe 2003js, 2006mk, 2006tk, and 
2007tt) the disagreement is larger than la, and explains the systematic negative offset between igpec and ^lc seen in 
Fig.Hni 
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